Simulation and test of the SLEGS TOF spectrometer at SSRF* 
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The Shanghai Laser Electron Gamma Source (SLEGS) is a powerful tool for exploring photonuclear physics, 
such as giant dipole resonance (GDR) and pygmy dipole resonance (PDR), which are the main mechanisms of 
collective nuclear motion. The goal of the SLEGS neutron time-of-flight (TOF) spectrometer is to measure GDR 
and specific nuclear structures in the energy region above the neutron threshold. The SLEGS TOF spectrometer 
was designed to hold 20 sets of EJ301 and LaBr3 detectors. Geant4 was used to simulate the efficiency of each 
detector and the entire spectrometer, which provides a reference for the selection of detectors and layout of 
the SLEGS TOF spectrometer. Under the events of "Pb, implementations of coincidence and time-of-flight 
technology for complex experiments are available; thus, ^; and neutron decay events can be separated. The 
performance of SLEGS TOF spectrometer was systematically evaluated using offline experiments, in which the 


time resolution reached approximately 0.9 ns. 
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I. INTRODUCTION 


The Shanghai Laser Electron Gamma Source (SLEGS), 
which is a high-intensity y-ray beamline based on laser 
Compton scattering (LCS) between 3.5 GeV, 200 mA rel- 
ativistic electron bunches in the storage ring of the Shang- 
hai Synchrotron Radiation Facility (SSRF) [1], and 10.64 um 
100 W CO; laser, was commissioned at the beginning of 2022 
[2. 3]. Owing to the innovative design of the interaction 
chamber [4] and dual collimation system [5, 6], SLEGS can 
operate in two modes to obtain tunable high-intensity Comp- 
ton spectrum of ^-rays from 0.66 MeV to 21.7 MeV with a 
full-spectrum flux of 10? - 10" photons/s. 

In addition, the CO; laser of the SLEGS can be operated 
in pulse mode with an adjustable laser interval and pulse 
width, for example, a 1 ms laser interval and 100 us pulse 
width at 10 W with a 1096 duty cycle. In the SSRF top-up 
and continuous-wave laser modes, the maximum beam flux is 
obtained, which is suitable for small cross-section measure- 
ments with some targets. More detailed information concern- 
ing the interaction chamber, dual collimation system, flux at- 
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tenuator, and commissioning results has been published else- 
where [4—7]. 


As a nuclear physics experimental platform, SLEGS aims 
at basic and applied research related to photonuclear reac- 
tions ranging from photon excitation to photon disintegration, 
such as nuclear resonance fluorescence, neutron production, 
and charged particle production [3, 8-10]. To achieve the 
above basic research goals, SLEGS has been equipped with 
four types of experimental spectrometers: gamma nuclear 
resonance fluorescence (NRF), neutron flat-efficiency detec- 
tor (FED) [11], TOF, and light charged-particle (LCP) spec- 
trometers, where both TOF and FED spectrometers can be 
used for (^, n) reaction measurements. (^, n) reactions have 
been measured using the ring-ratio method [12, 13] and direct 
neutron multiplicity sorting method [11, 14] using an FED 
spectrometer. However, FED spectrometers fail to measure 
neutron energy; therefore, it cannot acquire information on 
neutron energy. Similar to other TOF measurements [15-17], 
SLEGS TOF spectrometers can measure fast neutron energy, 
which is necessary for the study of isovector giant dipole res- 
onance (IVGDR) [18-21]. Because of the limitations of the 
construction and development of gamma sources, measure- 
ments of GDR based on photon sources have been restricted 
over the past 20 years [22]. 


The y-rays of SLEGS are currently available for experi- 
mental research on IVGDR. The SLEGS TOF spectrometer, 
which serves the same purpose as the ELIGANT [23, 24] at 
Extreme Light Infrastructure-Nuclear Physics (ELI-NP), was 
designed to measure the ^/-rays and neutron radiation emitted 
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Fig. 1. (Color online) A schematic layout of the SLEGS beamline. 


Table 1. Related parameters of the SLEGS TOF spectrometer. 


Detector Number Material Density (g/em®) Distance (cm) Diameter (inches) Thickness (inches) 


LaBr3 8-20 
EJ301 20 


LaBr3 
xylene 


5.10 
0.86 


30 3 4 
150 5 2 


during the decay of nuclear collective states. Therefore, the 
SLEGS TOF spectrometer has two main parts for y-ray and 
neutron detection. In addition, two mechanical supports were 
designed to hold 20xEJ301 and (8 — 20) xLaBrs detectors, 
as shown in Fig. 2. 

In this study, we focused on the simulation, calibration, and 
other tests of the SLEGS TOF spectrometer, which plays an 
essential role in photo-neutron measurements. 


Il. SIMULATION OF THE SLEGS TOF SPECTROMETER 


All simulations were based on Geant4 version 11.0.1 [25]. 
Various influencing factors were considered, and different in- 
teraction models were used for accurate simulations. For neu- 
tron interactions, QGSP_BERT_HP at neutrons < 20 MeV 
used high-precision neutron models and cross-sections to de- 
scribe elastic and inelastic scattering, capture, and fission. A 
G4NDL database was required for this physics list. More- 
over, RadioactiveDecay was activated. 

For the electromagnetic interactions, QGSP_BERT_HP 
was used as the default EM constructor. The emstan- 
dard_opt3, which is more suitable for the simulation of 
SLEGS TOF designed for any application, required a higher 
accuracy of electrons, hadrons, and ion tracking without a 
magnetic field. Therefore, emstandard_opt3 replaced the 
default EM constructor in the QGSP_BERT_HP. In this 
case, e*/e~ pair production was implemented using the 
BetheHeitler model with an LPM effect at high energies, 
and Compton scattering was implemented using the Klein- 
Nishina model. The photon-electric effect and Rayleigh scat- 
tering were both handled by the Livermore models. Different 
particle generators were constructed for different simulation 
purposes, such as monochromatic and isotropic gamma and 
neutron point sources. The geometries of the detectors and 


other parts were constructed consistently with the reality. The 
material properties of the detectors are listed in Table 1. 


A. TheSLEGS TOF spectrometer 


Fig. 2. (Color online) Layout of SLEGS TOF spectrometer. The 
LaBr; and EJ301 detectors(with two blue dotted ellipses) were 
placed at 30 cm and 150 cm from the target, respectively, all in the 
backward direction of the beam(with a red arrow). 


This section presents the construction of a complete 


SLEGS TOF detector array model, considering the factors 
that affect the measurement. (8 — 20) x LaBr3 and 20x EJ301 
were placed 30 cm and 150 cm from the target. The spec- 
trometer was placed close to the beam exit, far from the beam 
dump, and reserved for beam-spot monitoring, and an ex- 
ternal collimator was placed at the front of the experimental 
hutch. All detectors were placed in the backward hemisphere 
to reduce the background, as discussed in Sect. II D. Shield- 
ing, attenuators, and PMTs were also considered. 


B. Gamma detectors efficiency 


For GDR measurements with small cross-sections, the de- 
tection efficiency is an important consideration, which di- 
rectly determines the design of the experimental scheme and 
beam time. Gamma detectors are used to capture the y-rays 
emitted during the decay of collective states. In particular, 
when the SSRF is in the top-up mode and the laser is in con- 
tinuous output mode, LaBr3 must provide the starting time of 
TOF measurement by detecting de-excitation of y-rays due to 
neutron decay. 


Therefore, the gamma detectors must be fast and large 
enough to completely stop the y-rays produced by high- 
energy beams and close to the target for high absolute effi- 
ciencies and short time delays. A large-volume scintillator 
was selected because of its high efficiency, large size, and 
ease of maintenance. Considering its cost, LaBr3 detectors, 
which are used for beam monitoring and diagnosis, were cho- 
sen for gamma detection. Larger diameter detectors signifi- 
cantly increase the solid angle of coverage, while also imply- 
ing greater sensitivity to the beam background. Crystals with 
a diameter of 3 in are easily available. In [26], the charac- 
teristics of energy resolution, pulse shape, and time response 
of large volume LaBr; are discussed. The performance of 
LaBr3 with different thicknesses was simulated using attenu- 
ators and shielding. Fig. 3(a) shows the full energy peak in- 
trinsic efficiencies of single LaBr3 crystals of different sizes 
at y-ray energies below 22 MeV. 


The simulation results showed that a larger detector volume 
can provide higher intrinsic efficiencies. For a good balance 
between cost and y-ray efficiency, crystals with 3 in x 4 in 
LaBrs were used as the y-ray detectors in the SLEGS TOF 
spectrometer. Figure 3(b) shows the full energy peak abso- 
lute efficiencies of the 20 LaBrs arrays with respect to differ- 
ent y-ray energies. In this energy region (1 to 22 MeV), the 
full energy peak absolute efficiencies of the LaBrs array de- 
creased rapidly with increasing y energy. It can be seen that 
the full energy peak absolute efficiencies were less than 0.5% 
for the y-ray energy higher than 10 MeV. From the results 
of the peak-to-total ratio shown in Fig. 3(c), it can be seen 
that to obtain one 9 MeV ^-ray accumulated in the full en- 
ergy peak, the detectors must capture at least approximately 
30 y-rays on average. 
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Fig. 3. (Color online) (a) Full-energy peak intrinsic efficiencies of 
a single detector as a function of incident y-ray energies simulated 
for a LaBrs crystal with different sizes (see legend), (b) full-energy 
peak absolute efficiencies of 20x LaBrs with 3 in x4 in, and (c) peak 
total ratio of 20x LaBr3. Here, the monotonous energy point source 
generators were used, with y-ray energies from 1 MeV to 22 MeV. 


C. Neutron detectors efficiency 


Neutron detection and background control are important 
for nuclear astrophysics and technology applications [27—30]. 
In the IVGDR, decay neutrons may be also emitted, which 
can normally be measured using the TOF technique. The un- 
certainty in the distance and time resolution of the detector 
determines the accuracy of the reconstructed neutron energy. 
In general, a longer neutron flight distance is beneficial for 
reducing distance uncertainty. However, the solid angle cov- 
ered by the detector is inversely proportional to the square 


of the distance, resulting in lower absolute efficiency. There- 
fore, the time response and the efficiency of the detector are 
important reference factors. Organic liquid scintillation de- 
tectors are commonly used for the detection of fast neutrons. 
Owing to the abundance of hydrogen atoms in scintillators, 
the deposited energy of neutrons via recoiling protons mostly 
excites the electrons of scintillator molecules, and their de- 
excitation leads to the emission of optical photons. The clas- 
sic liquid scintillators BC501A, EJ301, and NE213 have been 
used in different arrays such as ENDA and ELINPGIANT- 
GN[3 1-33]. 
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Fig. 4. (Color online) (a) Simulated detector intrinsic efficiencies for 
an EJ301 scintillator detector of different sizes of 1 MeV to 6 MeV 
fast neutrons and (b) 20xEJ301 array absolute efficiencies in three 
positions: 100, 150 , and 200 cm from target. The mono-energetic 
isotropic neutron point source generators were used. 


However, EJ301 has been adopted as a substitute for 
BC501A, which exhibits almost the same performance [34, 
35] for neutron detection. Figure 4(a) shows the detector in- 
trinsic efficiencies simulated for EJ301 scintillator detectors 
of different sizes of 1 to 6 MeV fast neutrons. To increase 
the solid angle that the neutron detectors coverered, 20 EJ301 
detectors of 5 in in diameter were placed 150 cm from the 
target that suffered from an experimental hutch size restric- 
tion, which was irrevocable. A length of 2 in was consid- 
ered a good balance for efficiency, pulse shape discrimination 


(PSD), neutron flight distance uncertainty, and cost. 

Similarly, the neutron detection absolute efficiencies of 20 
EJ301 arrays placed 100, 150, and 200 cm from the target 
were also simulated, and the results are shown in Fig. 4(b). 
It should be emphasized that the variation at 4 MeV corre- 
sponded to the scattering cross-section of carbon [36]. The 
kinetic energy of a neutron is proportional to the square of its 
velocity, which can be determined by dividing its flight dis- 
tance by its flight time. Thus, the neutron energy resolution 
was determined by two factors: uncertainty of the neutron 
flight distance and uncertainty of the neutron flight time. For 
the former, the detector length and distance between the de- 
tector and target were determined together. For EJ301 with 
a thickness of 2 in, the uncertainties of the neutron flight 
paths were 2.43 46, 1.64 96, and 1.23 % at three different dis- 
tances of SI100, SI150, and SI200 cm, respectively. Ob- 
viously, longer flight paths resulted in more accurate neutron 
energy measurements, whereas the absolute efficiencies were 
reduced considerably, as shown in Fig. 4(b). The uncertainty 
of neutron flight time is discussed in Sect. III D. 


D. Background contribution 


Figure 5(a) shows the ^rays that hit the Pb target and pro- 
duced the aforementioned background. The most striking re- 
sult was that the emitted background did not scale with the 
target thickness. Moreover, particle trash multiple scatter- 
ing was evident. However, because the distance from the 
LaBr3 detectors was greater than 3 m, the absolute efficiency 
was extremely low; therefore, its influence was negligible. It 
is important to note that the mechanisms involved in high- 
energy y-rays were mainly Compton scattering and pair pro- 
duction, and the cross sections were several orders of mag- 
nitude higher than that of GDR. The simulation results re- 
vealed that when 9 MeV ^-rays interacted with a 1 mm lead 
target, the electron pair production process exceeded 70%, 
and the Compton scattering was close to 3096. Based on the 
simulation of the background components, the source of the 
background mainly came from ^-rays, whereas most of the 
charged particles, such as electrons and positrons, were ab- 
sorbed by the shieldings and attenuators of the detector. 

Figure 5(b) shows the angular distribution of the back- 
ground. The background emission was clearly forward; there- 
fore, the detectors were placed in the rear-angle direction. In 
particular, in the 90? direction, the minimum was reached ow- 
ing to the shadow effect and the self-absorption of the finite- 
size target. Figure 5(c) compares the background spectra sim- 
ulated in different directions for 9 MeV y-rays focused on a 
1mm Pb target, considering the hutch shielding wall and par- 
ticle trash. 


E. Pb target simulation 


Finally, the GDR of the Pb targets was used for the valida- 
tion and benchmarking of the code and the feasibility of co- 
incident measurement. The experimental data for GDR[22] 
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Fig. 5. (Color online) (a) Background energy spectra simulated for 20 x LaBr; detectors after emission of 9 MeV »^-rays into 0.1, 1, and 
5mm Pb targets, where the influence of experiment hutch shielding wall and particle trash were considered for some cases(see legend), (b) 
angular distribution of background simulated for 17 x LaBrs detectors placed 70 cm from target, 9 MeV ^-rays (10°) focused on a 1 mm Pb 
target, and (c) energy spectra of LaBrs simulated at different directions, with conditions shown in the panel. 


were fewer than those for calculations based on theoretical 
models [21, 37-41]. The reliability of the SLEGS TOF spec- 
trometer was verified using the GDR experimental data in 
208 

Pb. 
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Fig. 6. Partial decay scheme of ??5 Pp. 


To prove the feasibility of measuring the GDR using the 
SLEGS TOF spectrometer, a specific event generator in 
208Pb that encompassed seven channels of probabilities of 
decay[22, 32] was required: P1, P2, P3, and P4 cases were 
pure gamma decays, accounting for 2%, whereas the domi- 
nant decay was by emitting neutrons. In particular, the P6 
and P7 cases were suitable for measurement with SLEGS 
TOF spectrometer, where the simultaneously emitted ^-rays 
only required approximately 1 ns to travel through a flight dis- 
tance of 30 cm, whereas the neutrons travel through a flight 
distance of 150 cm to improve the accuracy of neutron en- 
ergy measurement. The time the EJ301 liquid scintillation 
detectors started determining neutron energy was the time at 
which the LaBr; detector placed at 30cm captured ^-rays 
used as. The seven different channel probabilities are sum- 


marized in Table 2 and are shown in Fig. 6. Based on the 
decay scheme, a generator comprising an isotropic neutron 
source and y source was coded. 


Table 2. Parameters of different decay channels[22, 32]. 
Case Particle Energy (MeV) Branch ratio Probability 


P1 ^y 9.00 93.46% 
yl 4.92 
P2 108 1.87% 
y6 6.39 
P3 q5 2.61 3.7496 296 
y3 5.81 
P4 44 0.58 0.9396 
y5 2.61 
P5 n3 1.60 35% 
n2 1.03 
P6 9 0.569062 ee - 
nl 0.70 
P7 78 0.32768 43% 
y9 0.56962 


The time spectra of the LaBrs detectors without any con- 
ditions are shown in Fig. 7(a). An obvious y-ray peak and 
three neutron peaks can be observed in the figure, where the 
y-ray peak consists of pure gamma decay and y-rays emitted 
concomitantly by neutron decay. For the three neutron peaks, 
the neutron decays corresponding to 1.6 , 1.03, and 0.7 MeV, 
were magnified and marked in the inset. The energy spectra 
when different time gates were used are shown in Fig. 7(b). 
The separated prompt y-ray peaks, delayed neutron events, 
and some background ^-rays can be clearly observed. The 
branching ratio was obtained from the relative intensity of the 
y-ray peak, and the corresponding energy level was obtained 
from the energy of the full energy peak or escaped peak under 
the appropriate time gate. In particular, the energy spectrum 
of LaBr3 shown in Fig. 8, with a 0-2 ns self-prompt time 
gate and y-ray peaks, is marked. 

Similarly, the time spectrum of EJ301 neutron detector 
without a time gate is shown in Fig. 9(a), whereas the en- 
ergy spectra with different time gates are shown in Fig. 9(b). 
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Fig. 7. (Color online) (a) Time-of-flight spectrum of LaBr3 without any gates. Three neutron peaks are magnified and their energies are 
marked in the inset. (b) Energy spectra of LaBr3 with some different time gates. The black lines represent the pure gamma decay and y- 
rays accompanied by neutron decay. The green line represents the background. The pink lines represent the neutron events and associated 


background. The resolution of the LaBr3 detectors was not considered. 
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Fig. 8. Energy spectrum of LaBr3 with a 0-2 ns self-prompt time 
gate condition. The resolution of LaBr3 was not considered. 


As expected, a clear separation of the neutrons decaying into 
different energy levels was achieved. 


F. Coincidence simulation results 


Using TOF technology for measuring the neutron energy, 
the energy resolution only reached an order of 10%. When 
the incident energy exceeded the double-neutron separation 
energy, more reaction channels were opened. It was impossi- 
ble to identify individual exclusive decays of GDR in heavy 
nuclei solely through neutron detection because of the limited 
energy resolution and dense-level schemes of residual nuclei. 
Therefore, the practicality of the coincidence technology im- 
plemented between different detectors was studied. The fold 
event distributions of the LaBr3 detectors, EJ301 detectors, 
and entire SLEGS TOF spectrometer are shown in Fig. 10. A 
fold event was a similar event recorded by multiple detectors. 


Most of the events were single. The multiple events detected 
by the LaBr3 detectors were significantly more than those 
detected by the EJ301 detectors because they were placed 
adjacent to each other and the generators involved multiple 
gamma particle emissions from GDR to low-lying states. In 
contrast, the EJ301 detectors were more dispersed because of 
limited budget. Multiple events were easily applied to select 
interesting and correlated events. The y — ^y coincidence en- 
ergy spectrum for LaBr3 is shown in Fig. 11(a) correspond- 
ing to multiple-step decays, such as the P2, P3, and P7 cases. 
In particular, an additional condition £1 -- F3 —9 MeV was 
applied to select two-step pure gamma decays, and the energy 
spectrum is shown in Fig. 11(b). The energy spectra of coin- 
cidence between ^-rays detected by LaBr3 and neutrons de- 
tected by EJ301 are shown in Fig. 12(a), which corresponded 
to the multistep decay of the P6 and P7 cases. These cases 
were the initial experiment of interest for the SLEGS TOF 
spectrometer. Because the start time signal was obtained by 
the LaBrs detectors, which was independent of the time signal 
from the beam, the complexity of the experiment was signifi- 
cantly reduced. In conclusion, itis feasible to select interested 
decay channels using the coincidence techniques. 


IIl. EXPERIMENT TEST RESULTS 
A. Data acquisition system 


The location of each component of SLEGS is shown in Fig. 
1. The SLEGS TOF was installed in the experimental hutch 
located between the laser and data acquisition hutches. In the 
experimental hutch, various spectrometers were used to con- 
duct the experiment using various methods. Some computers, 
which were connected to the laser system, beamline modula- 
tion system, and data acquisition plug-in via optical fiber con- 
trolled the beamline and acquired the spectrometer data in the 
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Fig. 9. (Color online) (a) Time-of-flight spectrum of EJ301 without a gate, and (b) energy spectra of EJ301 with different time gates. The 
black lines indicate the neutron events. The blue lines represent the pure gamma decay and y-rays accompanied by neutron decay, and the 
pink line represents the background. The resolution of the detectors was not considered. 
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Fig. 10. (Color online) Fold event distributions for LaBrs, EJ301, 
and entire SLEGS TOF spectrometer. 


data acquisition (DAQ) hutch. The walls of all hutches con- 
sisted of concrete, polyethylene, and steel. A particle trash 
that was designed to absorb neutrons and y-rays was built at 
the end of the experimental hutch. Therefore, an experimen- 
tal hutch is suitable for neutron-production experiments. 

The entire SLEGS TOF spectrometer data-acquisition sys- 
tem consisted of two parts: high-voltage control modules 
and DAQ. The multiple N1470 boards[42], which had four 
output channels, each capable of up to 8 kV (with a res- 
olution of +1V), provided high voltage for LaBr3 detec- 
tors. They were installed in a crate equipped with a com- 
pact UEP15-150 W[43] and a slot power supply, which pro- 
vided five direct current (DC) voltages: (+6 V 5A, +12 V 
3A, +24V 1.5A) The power supply of the EJ301 de- 
tectors was provided by double A7030DN boards[44] in- 
stalled in a SY5527LC universal multichannel power supply 
system[45]. Each A7030DN board had 12 channels, and the 
maximum output voltage and current reach 3 kV (with a res- 
olution of 1 V) and 1 mA (with a resolution of 1 uA), respec- 
tively. SY5527LC supported four slots with a maximum out- 


put power of 200 W (with 220 Vac Mains) and 350 W (with 
110 Vac Mains). In addition, the graphical control software 
and HiVoCS web tool performed remote operation via a host 
computer using TCP/IP communication protocol. 

The signals from each detector were sent to the DAQ via 3 m, 
a bayonet Nill-Concelman (BNC) cable, which had similar 
specifications to avoid attenuation of different channel sig- 
nals with a 50 Q impedance. The end of each cable was 
connected to a channel of digitizers via a BNC-micro coaxial 
connector (MCX) cable. Three CAEN V 1730s digitizers pro- 
vided 48 channels for 20x LaBr3 and 20x EJ301, which were 
installed in an NV8020 versa module europa (VME) crate. 
V1730s [46] was a 1-unit wide VME 6U module, housing 
a 16-channel 14-bit 500 MS/s FLASH ADC waveform digi- 
tizer with software selectable 2 Vpp or 0.5 Vpp input dynamic 
range on single-ended MCX coaxial connectors The DC off- 
set was adjustable in the +1 V (@2 Vpp) or +0.25 V (@0.5 
Vpp) range via a 16-bit digital-to-analog converter (DAC) 
in each channel. The module housed an optical link inter- 
face supporting a transfer rate of up to 80MB and provided 
a daisy chain capability. In addition to waveform recording 
firmware, CAEN provided four types of digital pulse pro- 
cessing(DPP) firmware: pulse shape discrimination (DPP- 
PSD), which is suitable for GDR measurement, pulse height 
analysis (DPP-PHA), zero length encoding (DPP-ZLEplus), 
and dynamic acquisition windows (DPP-DAW). Finally, a 40- 
m optical fiber was connected to a computer via an A3818 
peripheral component interconnect express (PCIe)[47]. The 
computer was installed with COMPASS readout software de- 
veloped by CAEN, which can set parameters, display wave- 
forms in real time, energy, timing, and PSD spectrum, and 
create root-format files. 


B. Energy calibration, resolution, and efficiency 


First, the existing 8 x LaBr3 energy and efficiency cali- 
brations were carried out using the Co and "Cs radiation 
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Fig. 11. (a) y — y coincidence energy spectrum of LaBrs, and (b) energy spectrum of LaBr3 with an additional gate of F1 + Ey2 =9 MeV. 


The resolution of the LaBr3 detectors was not considered. 
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Fig. 12. Energy spectra of coincidence between ^-rays detected by LaBr3 and neutrons detected by EJ301. In panel (a), coincidence y-ray 
energy spectrum is shown as a red line, corresponding to yg 327.68 keV and y9 569.62 keV, whereas coincidence neutron time spectrum is 
shown in panel (b) as a red line, corresponding to n» 1.03 MeV and n; 1.60 MeV. The resolutions of the LaBr3 and EJ301 detectors were not 


considered. 


sources, and the results are presented in Table 3. 

The full-energy peak efficiency comparison between the 
experimental data and the simulation is presented in Fig. 13, 
and it can be observed that they are almost consistent. Fig. 13 
shows the average experimental values of the full-energy peak 
efficiencies of the eight LaBr4 detectors. Based on the ex- 
periment with beam monitoring and diagnostics using LaBrs, 
the detectionr ranged from several hundred keV to 20 MeV. 
The detector measured synchrotron radiation up to approx- 
imately 300 MeV by reducing the high voltage. However, 
the detectors became blind when they faced high-flux gamma 
beams. In addition, because the LaBr3 detectors had excellent 
energy resolution and measurement range of energy, we mea- 
sured the direct y decay from GDR to the ground state and the 
low-energy ^-rays emitted from the GDR decay to low-lying 
states. 

The EJ301 liquid scintillation detector had poor energy res- 
olution, unlike the LaBr3 detector. Another difference was 


in the scintillation response. The response of EJ301 was 
nearly linear with the deposited energy for particles that re- 
sulted in electron recoil, such as gamma interaction with the 
scintillator[43]. However, serious nonlinearities can emerge 
in nuclear recoil events, such as neutron interaction with a 
scintillator [49-51]. Apparently, separate energy calibrations 
of y-rays and neutrons are necessary for EJ301 liquid scintil- 
lation detectors. In fact, EJ301 has been used to reconstruct 
neutron energy with TOF technology instead of directly deter- 
mining energy via pulse height analysis (PHA) and analog-to- 
digital converters (ADC). Therefore, only y-ray energy cali- 
bration of EJ301 was performed, and its time resolution is 
discussed in Sect. III D. 

As liquid scintillation is rich in hydrogen atoms, the Comp- 
ton scattering process was dominant in the interaction with 
y-rays. A full-energy peak was not observed in the energy 
spectrum; therefore, it is not feasible for the fitting method to 
perform energy calibration. The experimental and simulated 


Table 3. Parameters of eight LaBr3 resolution and efficiency calibration 


Model Threshold Cs 662 keV Co 1173 keV Co 1332 keV Calibration parameters 
AY. ) hyeshy FWHM Resolution Efficiency FWHM Resolution Efficiency FWHM Resolution Efficiency K B 

A18782-1 733 260.19 3.60 2.98% 45.76% 4.95 2.32% 35.45% 5.18 2.14% 31.40% 5.514 —4.479 
A18783-2 678 273.11 3.71 3.24% 45.19% 5.11 2.53% 34.46% 5.34 2.33% 30.98% 5.850 —7.686 
A18784-3 652 259.28 3.70 3.0996 45.62% 5.04 2.38% 32.97% | 5.24 2.19% 29.27% 5.598 —9.423 
A18785-4 744 269.65 3.62 3.12% 45.15% 4.88 2.37% 34.97% 5.06 2.17% 31.37% 5.743 —6.016 
A18786-5 638 245.84 3.65 2.87% 43.46% 4.90 2.18% 33.45% 5.06 1.9896 30.34% 5.247 — 6.013 
A18787-6 713 289.63 3.52 3.2196 44.17% 4.71 2.42% 34.55% 4.89 2.21% 30.57% 6.028 0.285 

A18788-7 696 250.03 3.50 3.02% 47.17% 4.46 2.23% 32.93% | 4.62 2.0396 29.4796 6.061 —40.900 
A18789-8 681 258.06 3.63 3.0196 47.42% 4.90 2.3196 32.01% 5.09 2.11% 28.15% 5.565 —9.062 
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Fig. 13. Comparison between experimental data using the Co and 
137Cs sources and simulation (see legend) for the full energy peak 
efficiency. The solid line connecting the simulation data is used to 
guide the eyes. 


energy spectra of "Cs are shown in Fig. 14(a). Energy de- 
positions from Compton back-scattered events as sharp edges 
were observed in the simulation plot, and events beyond that 
edge corresponded to multiple Compton scattering. In an ex- 
periment, the sharp Compton edge became smeared owing to 
the resolution effect. The simulation energy spectra of the 
137Cs and 6Co sources were broadened using Gaussian sam- 
pling to match the EJ301 detector resolution. The energy cal- 
ibration factors were obtained by comparing the experimental 
energy spectra with the simulated energy spectra after broad- 
ening. A comparison of the results is presented in Fig. 14(a). 
For the neutron detection efficiency of EJ301, a ?*?^Cf neu- 
tron source with an initial activity of 0.1531 uCi, which had 
a continuous neutron spectrum, was used. The *Cf source 
generally includes both alpha and spontaneous fission (SF) 
decay branches, with alpha decay comprising 96.91% of the 
activity and SF comprising the remaining 3.09% with neutron 
yield, which was approximately 4.4 x 10? neutrons per sec- 
ond per uCi[52]. Neutrons and y-rays were discriminated us- 
ing the PSD. The neutron detection efficiency was measured 
at different high PMT voltages, as shown in Fig. 14(c). For 
comparison, we simulated the neutron emission of the 252Cf 
source. The neutron spectrum was parameterized according 


to the following Maxwellian distribution[53]: 


N(E) = M Erorm-erm (1) 


where N(£) denotes the number of neutrons with energy E 
at a nuclear temperature of T' — 1.466 MeV. The simulated 
energy spectra obtained using this expression are shown in 
Fig. 14(b), and the efficiency was approximately 18%. This 
result was almost in agreement with the efficiency measured 
at the recommended high voltage of 1625 V. 


C. Pulse shape discrimination 


The EJ301 liquid scintillator detector was selected for neu- 
tron detection as part of the SLEGS TOF spectrometer ow- 
ing to its excellent PSD properties and fast timing perfor- 
mance. Because organic scintillators are sensitive to y-rays 
and neutrons, the pulse shape discrimination (PSD) capa- 
bility between neutrons and ^-rays was a critical feature 
that originated from three different components of scintilla- 
tion light, with mean decay times of 3.16 ns, 32.3ns, and 
270 ns[54, 55]. Different complex PSD methods have been 
reported, where the charge comparison (CC) method[57] 
based on the work of Jordanov and Knoll[56] and the zero- 
crossing method[58] are the most popular. In this section, 
the CC method, which is based on a comparison between 
the integrations of signals over two different intervals (entire 
pulse and only the tail) and the ratio of these two integrations 
as a separation parameter, was used to estimate EJ301 PSD 
performance. The following three critical parameters were 
written in firmware registers: PreGate, which determined the 
time before the pulse started to integrate, ShortGate, which 
determined the duration of the integration of the tail, and 
LongGate, which determined the total integration of the en- 
tire pulse. Therefore, PSD was introduced as a parameter for 
the implementation of the CC method, which is defined by 
the following equation: 


Qiong 


Qiong —Q short Corresponds to the integral of the tail, whereas 
Qiong and Qshort correspond to the integrals of the total pulse 
and Short Gate, respectively. Figure 15 shows the PSD 
spectrum of the *7Cf radioactive source measured using the 


PSD = (2) 
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Fig. 14. (Color online) (a) Energy spectra of "Cs source after calibration, initial spectrum, and Gaussian broadened spectrum based on 
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Fig. 15. (Color online) (a) 2D PSD spectrum of the 232 Cf radioactive source measured using the EJ301 detector, which was 100 cm from the 
source. (b) PSD histogram fitted using a double Gaussian function. The FWHMs and distance between y-ray and neutron peaks are marked. 


The FOM was equal to 0.6982. 


EJ301 detector. Pulses of neutrons generally have longer tails 
than those of y-rays, resulting in a larger PSD. Two distinct 
peaks were formed on both sides of the PSD, correspond- 
ing to y-rays and neutrons. Another crucial parameter was 
the figure-of-merit (FOM), which represents the resolution of 
neutrons and ^-rays, as defined by 


S 
FWHM, + FWHM, 


FOM = (3) 


where S is the separation between the peaks shown above, 
FW H M, is the full-width half-maximum of the spread of 
events classified as y-rays, and FW H M,, is the neutron 
peak. Some optimal setting parameters were necessary for 
high quality discrimination. Therefore, the effects of the 
Short Gate, Long Gate, Threshold, and high voltage were 
evaluated, and the FOMs calculated under different condi- 
tions are shown in Fig. 16. The high voltage and threshold 


mainly affected the detection efficiency rather than the FOM, 
whereas the Short Gate of 18ns and Long Gate of 200 ns 
were more appropriate for y/n discrimination. 


D. Timeresolution 


TOF technology was used to measure the fast neutron ener- 
gies. The time signal of the LaBr3 detectors or laser was used 
as the starting time and the time signal of the EJ301 detectors 
was used as the stopping time. Therefore, the time resolution 
of the detectors directly determined the accuracy of neutron 
energy measurements. 

For the time resolution of the LaBr; detectors, a pair of 
similar LaBr3 detectors performed the y-rays from the ??Cf 
source coincidentally. Additionally, a time resolution of 20 
x LaBr3 was studied in [23], which showed that single de- 
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Fig. 16. (Color online) FOMs calculated under different conditions. 
It should be noted that the integration gates started from a PreGate, 
and LSB was the least significant bit, which was converted to volt 
(1 LSB = (input dynamic range in Vpp) / 2%"). Here, 1 LSB = 
0.12mV. The solid lines connecting the experimental data were 
used to guide the eyes. 


tectors almost have similar time resolution. The coincidence 
time spectrum is shown in Fig. 17(a) The total FWHM 
was 1.28ns. According to equation AT. = ATi arsi? + 
AT apr 2^, assuming two detectors have similar performance, 
time resolution of a single LaBr3 was 0.90 ns. 


Similarly, a coincidence was observed between LaBrs and 
EJ301. To make the time resolution more reasonable, we 
sliced the long gate integration interval to 100 ch, and the 
experimental time resolution of each interval is shown in Fig. 
17(b) The uncertainty was taken as the root-mean-square 
(RMS). The average time resolutions for the 30 and 100 cm 
distances were 1.3 and 1.2 ns, respectively, and these contri- 
butions were from LaBr3 and EJ301. Considering that the 
time resolution of the LaBr; detector was 0.9 ns, the time 
resolutions of the EJ301 detector were 0.94 ns and 0.79 ns 
for distances of SI30 and SI100 cm, respectively. The re- 
sulting timing uncertainties were significantly smaller than 
the distance uncertainties of 2.43%, 1.64%, and 1.23% for 
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IV. CONCLUSION 


This study demonstrated that y-rays cover the vast majority 
of induced background, which cannot be absorbed by shield- 
ing. Therefore, a thinner target is preferable for reducing the 
induced background. Fortunately, regardless of the thickness 
of the target, the background-free signals of GDR decaying to 
the ground state and neutron signals of interest can be easily 
separated using appropriate time gate, PSD restriction condi- 
tions, and coincidence owing to the excellent time and energy 
resolution of detectors as well as y/n discrimination. 

In conclusion, the apparatus demonstrated exceptionally 
good performance in collecting and distinguishing ^-rays and 
neutrons of GDR decay. This will broaden our understanding 
of nuclear structures and lead to progress in the measurement 
of photon-neutron nuclear reactions. 
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